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Naval  Postgraduate  School 
Monterey,  California 


TVC  JET  VANE  THERMAL  MODELING  USING 
PARAMETRIC  SYSTEM  IDENTIFICATION 

I.  INTRODUCTION 

Jet  Vane  TVC  Systems 

Thrust  vector  control  (TVC)  systems  offer  means  of 
flight  vehicle  trajectory  control  that  are  virtually  inde¬ 
pendent  of  external  forces.  Such  a  capability  is  frequently 
required  for  taoticsl  missiles,  as  well  as  spacecraft  launch 
vehicles,  when  the  relative  flow  past  external  lifting  sur¬ 
faces  is  insufficient  to  generate  the  necessary  control 
forces.  This  oosaonly  occurs  during  low-speed  operations, 
such  ss  at  launch  or  during  hovering  flight.  High  angle  of 
attack  flight  say  also  lead  to  aituatlona  in  which  conven¬ 
tional  lifting  surfacee  are  inadequate.  In  addition,  there 
are  occasions  when  external  steering  devices  are  infeasible 
from  a  design  point  of  view,  such  os  for  tube-launched 
devices. 

Several  sethoda  of  TVC  hsvs  been  developed  end  applied 
to  operational  and  experimental  vahiclea.  These  include 
soveble  noxxlea,  internal  fluid  injection  (secondary  in¬ 
jection),  and  oechanical  jet  deflection  ayatems.  Jet  vsne 
systems  fall  in  ths  latter  category  and  they  tend  to  be 
favored  for  voluae-liaited  applicationa  requiring  relatively 
low  actuation  torques,  largs  thrust  deflection  angles,  and 
rapid  response.  Jet  vanes  aay  also  be  used  with  relative 


eas#  to  generate  roll  torques.  The  application  of  jet  vane 
TVC  dates  back  to  the  rockets  designed  by  Goddard,  and  has 
extended  to  the  Redstone,  Sergeant,  Taloa,  Pershing,  and 
Algo  II  and  III  motors  Cl, 23,  as  well  as  several  instal¬ 
lations  in  smaller  tactical  rockets. 

Of  course,  there  are  disadvantages  accompanying  the 
selection  of  jet  vanes  for  TVC  purposes.  These  include 
thrust  looses  on  the  order  of  3-5%  with  undeflected  vanes 
C23.  In  addition,  the  attainment  of  relatively  high  thrust 
deflection  angles  may  lead  to  axial  thrust  losses  of  the 
seme  order  of  magnitude  as  the  resulting  side  force.  How¬ 
ever,  the  chief  problem  associated  with  the  use  of  jet  vanes 
is  the  large  thermal  loading  that. they  experience  aa  they 
are  required  to  opsrats  in  hot,  high-speed,  particle-laden 
floes.  This  problsm  lesds  to  design  limitations  so  that  jet 
vanes  are  often  restricted  to  short-duration  uae  in  motors 
with  lov-teaperuture  non -seta lixed  propellants. 

The  aerodynamic  <eide-for-3e  producing)  character  1st  lea 
of  jet  vanes  say  be  calculated  with  fair  certainty  on  the 
baele  of  inviecld  flow  theory  with  suitable  corrections  for 
viscous  effects  (2,33.  On  the  other  hand,  difficulties  that 
stes  free  the  severity  of  the  jet -vane  thermal  environment 
have  led  to  design  practices  that  are  based  largely  upon 
past  experience  and  out-and-try  methods.  Over-design  is 
therefore  inevitable,  with  virtually  no  capacity  for  design 
optimization.  Zn  order  to  exploit  the  several  advantages  of 
jet  vans  TVC  aysteaa,  thsrefore,  it  has  become  necessary  to 


build  reliable  data  bases  and,  to  the  extent  possible,  at¬ 
tain  a  fundamental  understanding  of  the  heat  transfer  char¬ 
acteristics  of  such  systems. 

This  need  has  led  to  the  work  undertaken  at  the  Naval 
Postgraduate  School  (NPS)  in  support  of  a  larger  program  at 
the  Naval  Weapons  Center  (NWC).  Previous  investigations  at 
NPS  have  included  applications  of  computational  fluid  dyna¬ 
mics  (CFD)  [4,33  and  wind-tunnel  tests  using  infra-red  ther¬ 
mography  C63.  A  summary  of  the  results  of  these  studies, 
together  with  an  overview  of  relevant  previous  works,  are 
contained  in  Ref.  [73.  Work  in  the  area  of  CFD  has  con¬ 
tinued  (63  and  recent  reeulta  (S3  have  ehawn  that  dynamic 
simulation  methods  hold  promiaa  in  further  identifying  the 
dominant  fectore  affecting  the  thermal  characteristics  of 
jet  vanem.  Thie  latter  area  of  study  has  advanced  to  the 
application  of  parametric  eyatem  identif iomtion,  and  it  is 
the  reeulta  of  these  efforts  that  are  the  main  subject  of 
this  report. 


In  the  design  of  a  jet  vane  eyatem,  the  integrity  of 
the  vane*  themselves  must  be  guaranteed  over  the  specified 
work  cycle.  Although  this  is  a  serious  challenge,  designers 
aust  also  consider  the  behavior  of  the  vanes  and  supporting 


thermal  stresses  may  develop  due  to  the  proximity  of  a  rela¬ 
tively  cool  supporting  structure.  These  and  other  design 
aspects  can  only  be  addressed  with  precision  if  there  is  a 
good  understanding  of  the  convective  heat  transfer  process 
that  gives  rise  to  the  energy  transfer  from  the  flowing 
gases  to  the  the  vane.  Put  another  way,  the  application  of 
complex  computer  codes  for  thermal  conduction  in  the  vane 
and  supporting  structure  can  only  follow  the  specification 
of  the  convective  boundary  oonditiona. 

With  respect  to  these  conditions,  the  problem  is  even 
more  complex.  The  vane  is  immersed  in  a  flow  field  that  is, 
if  generally  described,  compressible,  turbulent,  multi-com¬ 
ponent  (and  possibly  multi -phase),  three-dimensional,  and 
unsteady,  with  variable  properties  and  nonlinear  and  time- 
variant  boundary  conditions.  Even  if  taken  one  at  a  time, 
the**  complexities  present  problems  thet  are  beyond  the 
state-of-the-art  for  exeet  solution.  In  addition,  the  over¬ 
all  flow  field  will  contain  intersecting  and  impinging  shock 
waves  that  give  rise  to  discontinuous  evsntm  and  further 
complication  of  tha  boundary  oonditiona.  The  presence  of 
various  protuber -sncsa  only  serves  to  exacerbate  thee* 
difficulties. 

To  further  define  the  problem,  it  mmy  be  of  use  to 
consider  the  levels  of  heat  transfer  thet  might  be  expected 
from  the  point  of  view  of  "simple*  convection  from  a  super¬ 
sonic  flow  to  e  cooled  wall.  The  convective  heat  transfer 
coefficient  <h>  may  be  described  in  terms  of  the  nondimen - 

t.j  < 


aional  Stanton  number  (St)  aa  follows: 


h  =  (pVc  )  St 
P 

where  p,  V,  and  are  the  density,  velocity,  and  specific 
heat  at  constant  pressure  of  the  flowing  gas.  On  the  as¬ 
sumption  that  the  gas  behaves  ideally,  this  expression  may 
be  written 

h  ■  St  ‘iTT  T  8  ImV> 

where  P  and  T  are  local  gas  pressure  and  temperature  and  I 
is  the  specific  impulse.  Consider,  for  instance,  the  condi¬ 
tions  at  the  exit  of  a  rocket  nozzle  flowing  at  Mach  3. 

With  a  ratio  of  specific  heate  assumed  to  be  k  »  1.2,  we 
have,  approximately: 

P 

h  ■  St  (0.25  ^  g  I  ) 

To  *P 

where  the  subscript  <q>  refers  to  stagnation  conditions.  If 
a  typical  case  is  taken  to  be  given  by  Pq  *  100  bar,  Tq  » 

2 St?0  and  I  _  "  250  a,  then,  in  round  nuabers,  the  heat 

.  v  v  -  .  ttp 

S. 

transfer  coefficient,  in  units  of  (kW/a*K>,  is  h  =>  2000  St. 

In  rocket  engine  nozzles,  a  typical  value  for  the 
Stanton  number  is  about  0.002  1101,  Thus,  vlth  the  result 
above,  the  expected  value  of  the  heat  transfer  coefficient 
fa  on  the  order  of  4  kW/»* K,  Such  a  value  should  not  be 
taken  as  conservative,  alnoe  actual  motor  conditions  may  be 
sore  severe  and  local  regions  (e.g.  stagnation  points  or 
regions  of  flow  reversal)  may  exparlence  such  higher  values. 
By  way  of  comparison,  values  of  h  such  as  these  are  much 


larger  than  those  usually  considered  for  convective  heat 
transfer  and,  in  fact,  they  are  more  typical  oi  those  real¬ 
ized  in  phase-change  processes  C113. 

II.  MODELING  AND  SIMULATION 

Background 

A  general  goal  of  jet  vane  heat  transfer  studies  is  to 
develop  a  capability  by  which  the  transient  thermal  behavior 
of  the  vane  may  be  predicted  with  confidence.  An  under¬ 
standing  of  the  energy  transport  proceasea  occurring  at  the 
boundaries  of  the  vane  is,  of  course,  essential  to  the 
achievement  of  this  goal  but,  because  of  the  complexity  of 
the  thereof luid  environment  in  whioh  the  vane  operates,  CFD 
methods  cannot  be  expected  to  yield  a  complete  and  compre¬ 
hensive  predictive  capability.  In  addition,  the  results  of 
CFD  modeling  must  invevitably  be  simplified  if  they  are  to 
be  adapted  to  vane  design  purposes  and  the  prediction  of 
system  performance.  In  order  to  complement  the  CFD  studies, 
therefore,  a  program  has  been  initiated  to  approach  the 
problem  fros  a  lees-general  but  more-practical  point  of 
vies.  Underlying  this  approach  is  an  acceptance  of  the  fact 
that  details  of  the  energy  transport  proceasea  will  never  be 
fully  defined.  On  the  other  hand,  there  is  also  the  recog¬ 
nition  that  the  cumulative  or  'lumped*  effect  of  all  the 
various  complications  mentioned  above  is  to  transfer  energy 
to  and  fros  the  vane.  These  boundary  processes  lead  to 
further  energy  flow  within  the  vane  that  ultimately  gives 


rise  to  a  vans  temperature  distribution  that  varies  with 
time. 

This  transient  behavior  should  reflect  the  nature  of 
the  boundary  heat  transfer  processes  that  drive  it,  as  veil 
as  the  thermal  impact  of  various  vane  design  parameters 
(thermal  conductivity,  heat  capacity,  density,  geometry, 
etc.).  The  transient  response  is  the  'signature”  of  the 
combined  effects  of  the  thermal  environment  and  the  vane 
configuration.  This  premise  is  the  basis  for  justification 
of  the  modeling  and  simulation  (MIS)  study  described  in 
subsequent  sections  of  this  report. 

The  goal  of  the  M&S  etudy  has  been  to  develop  a  dynamic 
model  of  vane  heating  and  cooling.  If  such  a  model  can  be 
constructed*  even  though  it  is  approximate,  then  significant 
benefits  will  accrue.  For  instance,  measured  vane  tempera¬ 
ture  histories  can  be  used  for  deductive  purposes  ev.ch  as 
parametric  system  identification  (PSD  whereby  the  model  la 
used  to  deduce  what  would  have  had  to  be  true  in  order  for 
the  obssrved  temperature  history  to  have  occurred.  Having 
estimated  the  values  of  appropriate  parameters,  local  tem¬ 
perature#  can  be  deduced  at  points  in  the  vane  structure 
where  conditione  (accessibility,  sensor  eurvival,  etc.  > 
preclude  measurement.  In  addition,  a  trustworthy  predictive 
model  would i  (1)  reduce  the  occaslone  in  which  testing  is 
required  to  verify  the  effects  of  design  changes,  (2)  permit 
the  extension  of  sub-scale  test  data  to  predict  full  scale 
jet  vane  performance,  and  (3)  indicate  the  design  directions 


moat  likely  to  lead  to  optimum  system  performance. 

Approach 

The  HfcS  approach  may  be  thought  of  as  a  vast  simplifi¬ 
cation  of  a  numerical  model  that,  were  it  possible  to  con¬ 
struct,  would  lead  to  a  successful  predictive  capability. 
Whereas  such  a  comprehensive  model  would  treat  the  flow 
environment  of  the  vane  and  the  vane  itself  in  fine  numeri¬ 
cal  detcil,  the  HfcS  method  assumes  that  sufficient  accuracy 
may  be  obtained  if  the  flow  and  the  vane  are  made  up  of  a 
relatively  few  thermal  parts.  An  important  product  of  the 
HfcS  approach  is  the  prediction  of  the  time  history  of  the 
vane  temperature  -“in  fact,  this  is  what  provides  the  vital 
clues  that  are  used  to  estimate  the  values  of  the  model 
parameters. 

The  system  identif ieation  procedures  used  in  this  study 
are  those  faaillsr  to  the  discipline  of  automatic  controls. 
If  the  governing  sathesstiosl  relationships  can  be  cast  in  a 
transfer -function  or  other  eyetesstic  form  (not  neceeearily 
linear),  then  experimental  observations  can  be  related  to 
the  various  parameters  of  the  model  by  a,, «tea  identification 
methods.  In  order  for  thia  procedure  to  be  of  value,  how¬ 
ever,  the  afore -mentioned  parameters  must  be  related  to  the 
physical  quantities  affecting  the  thereal  transport  process. 
The  establishment  of  these  relationships  is  described  in  the 
nest  section  of  thie  report. 


III.  MODEL  DISCRETIZATION 


The  model  as  presently  configured  considers  only 
thermal  energy  input  processes:  forced  convection  at  the 
vane  surface  and  stagnation-point  heat  transfer  at  the 
leading  edge.  These  processes  are  driven  by  the  flow  stag¬ 
nation  temperature  wnich  is,  in  turn,  derived  from  the 
measured  thrust  levels  by  means  of  the  rocket -motor  bal¬ 
listic  characteristics.  The  rise  and  fall  of  the  stagnation 
temperature  ia  assumed  to  follow  that  of  the  thrust  without 
significant  dynamic  offset.  In  other  words,  the  generic 
product  of  the  model  is  a  aeries  of  functions  that  give 
local  vane  temperatures  in  response  to  thrust.  (If  these 
functions  were  linear,  they  would  have  direct  transfer- 
function  counterparts.  > 

Several  vane  discretization  schemes  hsve  been  tested, 
and  it  has  been  found  that  four  'lumps*  are  sufficient  to 
indicate  the  general  thermal  behavior  of  the  vane:  vane  tip, 
vane  body,  shaft,  and  mount.  This  conf iguration  leads  to  a 
simplicity  that  is  s  necessary  feature  of  the  H4S  approach, 
end  resulte  to  be  described  later  justify  the  use  of  euch 
simple  models  ss  design  tools. 

Estimation  p£  Geometrical  and  Thermal  Properties 

As  a  preliminary  step  in  the  atudy,  the  actual  vane 
design  of  interest  to  the  Naval  Weapons  Center  was  confi¬ 
gured  as  s  collection  of  lumps  with  geometries  suitable  for 
estimation  of  thermal  conduction  properties.  Figure  1  il¬ 
lustrates  the  vane  geometry  that  has  been  modeled  (the  vane 


la  constructed  of  1054Cu/W,  and  the  mount  material  is  steel). 

From  this  design,  the  fictitious  thermal  vane  has  been 
hypothesized  as  consisting  of  three  rectangular  solids:  tip, 
fin,  and  shaft.  This  discretized  vane  is  shown  in  Fig.  2. 
The  tip  of  the  van©  has  been  separately  identified  in  order 
to  account  for  the  stagnation  properties  of  the  thermal  con¬ 
vection  near  the  vane  leading  edge.  This  portion  haa  been 
arbitrarily  sized  so  as  to  have  a  chord  length  of  10X  of  the 
total  vane  chord  of  3.75  inches.  With  this  length  fixed, 
the  discretized  tip  is  shaped  bo  as  to  have  the  same  lateral 
area  (chord  x  thickness)  as  that  of  thrr  actual  tapered  and 
rounded  tip. 

The  remainder  of  the  discretized  fin  is  assumed  to  be 
subjected  to  thermal  convection  of  s  turbulent  boundary 
layer  type.  It  is  lumped  into  s  rectangular  solid  of  thick¬ 
ness  equal  to  the  average  value  of  the  tapered  fin  (with  tip 
removed)  end  spmn  equal  to  the  mean  span  of  the  actual  fin. 
The  remaining  length  dimension  is  set  by  equating  the  volume 
of  the  rectangular  solid  to  that  of  the  fin  (lees  the  tip 
portion).  The  van#  shaft  is  similarly  "molded6  into  a  rec¬ 
tangular  solid,  and  this  element  ie  ameumed  to  be  subjected 
only  to  conduction  (and  radiation,  were  it  included)  heat 
transfer. 

The  vane  mount  (not  ehovn  in  Flge.  1  12)  ia  a  rela¬ 
tively  massive  and  complex  structure.  In  order  to  pursue 
the  basic  feasibili  ty  of  the  method,  no  attempt  vaa  mad  ->  to 
model  the  thermal  resistance  of  this  component  in  detail. 
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instead,  representative  seta  of  values  (conduction  lengtn 
and  cross-section)  for  the  main  mount  components  were  esti¬ 
mated  using  a  rough  scaling  procedure.  The  appropriate 
thermal  resistance  of  the  mount  was  then  computed  from  an 
analysis  of  the  analogous  electrical  circuit.  In  prelimi¬ 
nary  studies  the  mount  thermal  resistance  was  considered  to 
be  an  •adjustable  parameter"  for  use  in  seeking  agreement 
**itv  the  eveilwtle  experimental  evidence.  (Determination  of 
the  correct  value  of  this  quantity  is  one  of  the  goals  of 
the  ssysiem  identification  procedure  described  later. > 

Each  of  the  c  mpo. tents  described  above  was  assigned  a 
thermal  node  located  at  it's  mass-canter  which,  in  turn,  is 
the  assumed  loc^ti^n  of  the  energy  storage  associated  with 
the  entire  maaa.  Ttv^  capacity  to  store  energy  vaa  calcu¬ 
lated  in  the  usual  way  for  each  lump  For  the  mount,  the 
thermal  cape -pity  was  eeeuued  to  be  infinite  --  thus  the 
temperature  of  this  node  remained  constant  at  ambient 
(ground)  temperature  during  the  simulation. 

It  should  be  emphasised  that  at  thia  point  in  the  atudy 
the?  discretization  rational#  ie  quite  arbitrary  --  che  goal 
ha*  boon  to  obtain  adequate  agreement  «>*. -h  teat  rerulte 
using  ttUD  minimum  number  of  thermal  components.  Refinement 
of  the  model  can  easily  include  tne  division  of  the  vane 
atructure  into  a  larger  number  of  emaller  "lumps."  Altnough 
thia  would  permit  the  estimation  of  temperatures  at  more 
nodes  and  intarfacea,  the  accuracy  of  these  estimates  would 
be  no  greater  than  that  associated  witn  the  "minimum  node" 


model  described  here.  Such  added  complication  of  the  vane 
nodal  distribution  would  only  be  justifiable  in  conjunction 
with  a  more-detailed  description  of  the  thermal  boundary 
conditions,  and  this  would  add  even  more  uncertainty  to  the 
model.  Results  described  later  indicate  that  such 
complication  is  not  necessary. 

Development  of  the  Governing  Equations 

With  the  thermal  network  and  associated  properties  thus 
defined,  the  governing  equations  were  formulated  by  means  of 
an  energy  balance  at  each  node.  In  general,  this  balance 
reads  "rate  of  heat  flow  in  3  rate  of  heat  flow  out  ♦  rate 
of  energy  storage. *  For  the  fin  node  (node  2  in  Fig. 3),  for 


instance} 


M 


T  T  -  T  T  -  1 
*2  i  1 2  o  J2 \ 

RF2  Ri2  R23 


a  .  C2sT2 


(  l  > 


Here  the  symbol  R  is  used  to  denote  thermal  resistance  -- 

1/hA  for  convection  to  the  fin  and  L/kA  for  conduction 

between  the  nodes  --  and  C  is  thermal  capacitance  (mass  x 

heat  capacity).  The  letter  a  is,  as  usual,  the  symbol  for 

the  Laplace  variable.  Calculation  of  the  recovery  temper- 

ature  T^  ie  discussed  in  the  next  section. 

The  temperature  at  the  node  may  be  expressed  explicitly 

in  tarns  of  the  surrounding  nodal  temperatures  as  follows i 

R 


T  T  T 

_ *.S.3 _  <152  .  Ll  .13 

1  *  t2s  Rp2  Ri2  R23 


) 


(2) 


where  the  nodal  resistance  Rn2  is  given  bys 
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Figure  3.  Nodal  configuration  and  table  of  estimated 
values. 


16 


and  the  nodal  time  constant  is  given  by 
a2  *  Rn2  C2 

Nodes  of  the  model  that  do  not  provide  an  energy  storage 
function  are  also  absent  the  associated  time  constant.  Thus 
the  temperature  at  the  internal  node  shown  at  the  interface 
between  the  tip  and  fin  of  the  present  model  (node  i  in  Fig. 
3)  is  expressed  in  terras  of  the  temperature  of  the  surround¬ 
ing  nodes  with  only  thermal  resistances  as  parameters. 

A  heat  balance  at  each  of  the  remaining  nodes  yields 
relationships  similar  to  those  given  above*  so  sufficient 
information  is  available  to  determine  each  nodal  temperature 
as  a  function  of  time.  Initially*  the  aystem  is  assumed  to 
be  at  thermal  equilibrium  with  the  environment  so  that  the 
nodal  temperatures  are  all  equal  to  the  ambient  value.  With 
the  firing  of  the  rocket  motor  this  equilibrium  is  disturbed 
and  the  two  recovery  temperatures  provide  inputs  to  the 
ensuing  heat  transfer  process.  In  preliminary  studies  the 
Digital  Simulation  Language  <DSL)  hae  been  used  to  automati¬ 
cally  accomplish  the  repeated  integration  steps  required  to 
determine  the  thermal  response  of  those  nodes  possessing 
storage  capacities. 

The  ualoulation  method*  together  with  preliminary 
results*  are  described  in  later  aeationa.  Before  this* 
howe^. .  it  is  necessary  to  estimate  the  driving  functions 
V »  entire  proaesa  -*  the  thermal  inputs  at  the  tip  and 
fin  aurfacea.  This  ie  the  subject  of  the  next  chapter. 


IV.  ESTIMATES  OF  SURFACE  HEAT  TRANSFER  RATES 


As  described  in  the  previous  section,  the  forcing  func¬ 
tions  of  particular  interest  are  those  due  to  heating  at  the 
vane  tip  and  along  the  downstream  surface  of  the  vane.  (The 
cooling  effects  of  radiation  and  ablation,  both  clearly  im¬ 
portant  in  the  actual  situation,  are  omitted  from  discussion 
for  the  time  being. )  Most  of  the  heat  transfer  computations 
outlined  below  are  derived  from  the  analytical  methods  sug¬ 
gested  in  the  works  edited  by  C. C.  Lin  £121  and  in  the  AGARD 
Monograph  authored  by  Ziebland  and  Parkinson  £10 

A  particular  feature  of  these  high  speed  flows  is  the 
large  difference*  in  temperature  that  the  gas  experiences  in 
decelerating  near  the  body  surface.  In  auch  cases  it  is 
necessary  to  account  for  the  temperature  dependency  of  gas 
properties.  A  maet-ueeful  simplification  for  this  purpose 
ie  that  proposed  by  Oieeeler  and  others  £12,  p.3043  in  which 
the  Prendtl  number  end  specific  heat  are  considered  to  be 
constant  inasmuch  as  their  variations  with  temperature  are 
of  a  lower  order  of  magnitude  than  those  of  the  other  gaa 
properties  (viscosity  and  thermal  conductivity).  Thus,  with 
this  assumption,  the  quantity  Pr/cp  «  constant  and  a  separ¬ 
ata  estimate  of  the  gas  viscosity  leads  to  the  thermal 
conductivity  for  a  given  Prandtl  number  and  mpeaific  heat. 

For  the  rough  calculations  used  in  this  study,  the 
Prandtl  number  has  been  estimated  using  the  Eucken  formula 
£2,  p. 139J,  Pr  »  4k/<9k-5),  where  k  ia  the  ratio  of  specific 
heats.  With  the  gas  constant,  Rg,  the  constant -pressure 


R  k/(k-i). 

g 


specific  heat  la  given  by  c^  =  R  k/<k-l). 

In  this  work,  the  standard  Sutherland-type  formula  has 
been  adopted  for  determining  viscosity: 

-  if-)" 

^ref  Tref 

Following  the  recommendation  of  C 2,  p.9],  a  value  of  n  =  0.7 


has  been  used.  In  addition,  the  reference  viscosity  has 


been  taken  to  be  p 


4.0x10  5  N-a/ra2  at  T  »  1000  K. 

ref 


In  the  analyaia  of  stagnation  point  heating,  the  solu¬ 
tion  to  the  boundary  layer  equation  requires  an  estimate  of 
the  local  fluid  acceleration  in  the  vicinity  of  that  point. 
In  the  case  of  supersonic  flow,  this  can  be  approximated 
upon  the  assumption  that  Newtonian  flow  prevails  between  the 
■  ■bow' wave'-. 

bow  Wave  and  the  body  and,  (see  van  Driest  Cl,  p* 3661), 

6  V  <du  /dx>  <U/t>>  lacp  /P  XT  /T  >  l°* 5 

m  X*Q  •  ay  Oa  » 

Hsire  the  subscript  <•>  refers  to  freeatreas  oonditi on«,  (o) 


denotes  stagnation  conditions,  and  (y)  is  for  conditions 
downstreas  of  the  normal  shock.  The  pressure  and  tempera¬ 
ture  ratios  in  ths  above  expressions  are  known  functions  of 
the  freestreaa  Haoh  number. 

Nith  8  thus  defined,  the  Stanton  number  may  be  computed 
from  the  following i 


St  *  0.57  taD/U>°‘SPr"0*6  Re~°'S 

a  o  a 
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Turbulent  boundary  Laver  Convection 

In  high  speed  compressible  flows,  the  analysis  of  ther¬ 
mal  pracesaes  is  complicated  by  the  fact  that  considerable 
compression  and  viscous  dissipation  follow  from  the  deceler¬ 
ations  occuring  in  the  boundary  layer.  This  leads  to  tem¬ 
peratures  within  the  boundary  Layer  that  are  in  excess  of 
that  of  the  freeatream,  and  the  driving  temperature  for  heat 
transfer  is  the  so-called  recovery  temperature,  Tn.  Thus  a 
recovery  factor  ia  defined  aa  follows: 


Fortunately,  it  has  been  found  by  many  investigators  that 
the  recovery  factor  may  be  related  to  the  Prandtl  number  in 
a  simple  way  that  is  adequate  for  meet  purposes.  For  a  tur¬ 
bulent  boundary  layer  (assumed  here),  the  relationship  is 


r  *  Pr 


1/3 


A  more-difficult  problem  arises  from  the  dependency  of 
gee  properties  upon  temperature  since  in  moat  aaaee  the 
appropriate  reference  temperature  for  thia  calculation 
depends  upon  the  well  temperature.  This,  in  turn,  depends 
Upon  the  reference  temperature.  The  reference  temperature 
is  typically  defined  ae  (23 

T»f  “  °-s  T«n  *  °'2e  T.  *  °'22  tr 

The  present  calculations  are  based  upon  the  assumption  that 
the  eel 1  temperature  ie  adequately  represented  for  this  pur¬ 
pose  by  the  seen  of  the  recovery  end  ambient  temperatures. 
(An  evaluation  of  the  validity  of  this  assumption  requires 
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an  iterative  time-dependent  calculation.  In  the  system 
identification  work  to  be  described,  however,  the  convective 
resistance  is  treated  as  a  parameter  to  be  identified.  Thus 
the  experimental  data  lead  to  the  deduction  of  the  "effec¬ 
tive1*  film  coefficient  and  the  detailed  analyses  described 
in  this  section  may  not  be  necessary  for  design  purposes. ) 

Given  that  the  reference  temperature  for  gas  properties 
is  adequately  described,  the  gas  viscosity  and  thermal  con¬ 
ductivity  may  be  estimated  aa  above.  The  Stanton  number  for 


turbulent  compressible  flow  is  then  given  by  C21: 


St  »  G.029G  Pr'0,67  Re“0*2 


From  the  expressions  given  above  for  the  Stanton  num¬ 
bers,  Huseelt  numbers  and  thermal  resistances  can  be  calcu¬ 
lated  for  the  model  components  affected  by  stagnation  and 
boundary  layer  heat  tranefer  processes.  It  should  be  noted 
that  the  thermal  resistances  are  scale-dependents  stagnation 

-15 

point  thermal  resistance  decreases  as  (scale)  .  ’ ,  and  turbu¬ 
lent  boundary  layer  thermal  resistance  is  propportionsl  to 
'  -i.  2 

(scale)  ,  The  simulations  conducted  in  this  study  have 
been  of  a  1 /4-scale  van*  in  order  to  provide  a  comparison 
with  available  HVC  test  data. 


The  information  necessary  to  provide  values  for  the 
thermal  resistance  includes  the  freestreso  characteristics 
at  the  location  of  the  Jet  vane  *•-  Heoh  number,  stagnation 
temperature,  and  stagnation  pressure.  The  concept  of  the 


present  model  ie  that  theae  quantities  are  determined  from 


z\ 


the  ballistic  characteristics  of  the  rocket  motor.  Accord¬ 
ingly,  inputs  required  for  the  simulation  include  the  motor 
chamber  pressure,  thrust,  and  characteristic  velocity,  and 
the  discharge  coefficient  and  pressure  ratio  of  the  nozzle. 
In  addition,  necessary  propellant  gas  properties  include  the 
ratio  of  specific  heats  and  the  gaa  molecular  weight  or  gas 
constant.  Experimental  values  are  preferred,  to  the  extent 
that  they  are  available,  but  a  considerable  body  of  theory 
exists  if  analytical  estimates  are  neceeaary. 

From  these  quantities,  the  stagnation  temperature  may 
be  calculated  as  follows i 

T  -  CC  .  r  c*)2/R 
a  d  g 

where : 

C<1  -  di.alu.rg.  ao.Mtel.nt  ..930 

r  •  ik<2/k.l)lk*1>/<k"1>l°‘s  (0.65) 

a*  »  characteristic  velocity  (1512  »/e> 

K  «  gas  constant  <318.5  »2/e2-K) 

k  «  ratio  of  specific  heats  < 1. 213 

the  numbers  given  in  parentheses  are  these  presently  in  use, 
and  lead  to  a  stagnation  temperature  of  2850  K. 

From  the  motor  nozzle  pressure  ratio  Mach  number  at  the 
exit  ( assumed  to  be  that  at  the  vane)  la  computed  from 

2  ( is  —  i  )  /  if 

M  *  C2y<k-l)l  C(P  /P  )  *  -  13 

S  Q«  m 

and  selecting  a  nozzle  pressure  ratio  of  p0-,/?<>  *  186,  this 
expression  gives  «  3. 75.  It  is  important  to  note  that 
the  motor  chamber  pressure  must  also  be  provided  because  the 
heat  transfer  calculations  for  the  vane  require  the  density 


and  hence  the  pressure  in  the  freeatream  at  the  vane  loca¬ 
tion  --  in  the  calculations  pi'eeented  here,  the  value  of 

P  »  IS. 76  MPa  has  been  used. 

Qm 

For  ease  of  reference,  the  following  table  of  values 

describes  the  heat  transfer  quantities  obtained  from  the 

inputs  and  computed  quantities  described  above. 

Stagnation  Turbulent  Boundary 
Point  Layer 


Stanton  number 
Mussel t  number 
Fils  coefficient, 
Thermal  resistance,  K/W 


а.  53x10 
43.  1 

б.  46xi03 


5.  69 


3.  01x10 


436 

2. 25xl03 
0.  634 


V.  SIMULATION  RESULTS 

The  model  and  governing  equation#,  described  above, 
have  been  coded  in  the  Digital  Simulation  Language  (DSL), 
This  tBN  so it were  product  is  a  FORTRAN -based  language  with 
the  particular  advantage  that  a  vide  variety  of  integration 
subroutines  are  imbedded  and  directly  available  to  the  utter. 
Such  integrations  are  necessary  at  each  of  the  energy  stor¬ 
age  node#,  and  for  this  purpose  a  5-th  order  Runge-Kutta 
method  Incorporating  a  variable  stop  size  was  selected.  In 
addition  to  the  vane  conduction  model  described  above,  the 
code  incorporetee  energy  flow  calculations  in  order  to  veri¬ 
fy  that  the  model  does  indeed  account  for  the  disposition  of 
aifl  heat  transferred  to  or  froa  the  vane,  me  well  ae  that 


ZJ 


stored  at  the  various  nodes. 

The  simulation  code  has  been  written  so  that  the  dri¬ 
ving  input  is  the  thrust  of  the  rocket  motor,  as  provided  by 
test  data.  In  the  case  reported  here,  che  point -by-point 
thrust  vs.  time  measurements  are  approximated  by  a  ramp  in¬ 
put  from  0  to  2335  N  <525  lh)  in  a  period  of  0.5s,  held  at 
this  ievel  for  2.5s,  and  ramped  back  to  zero  thrust  in  a 
further  period  of  Q.  5a.  Although  the  actual  thrust  data 
points  could  be  used,  the  rarap-up/rarap-dovn  closely  approxi¬ 
mate*  the  thrust  schedule  and  provides  a  useful  input  for 
evaluating  the  aensitivity  of  the  thermal  response  of  the 
vane  to  varioua  parameters  of  the  model. 

An  important  aspect  of  the  physical  ©vent  of  rocket 
firing  is  that  the  convective  resistance*  to  vane  heat 
transfer  tat  the  stagnation  point  and  in  the  boundary  layer) 
are  coupled  to  the  preaence  or  absence  of  ilov  past  the 
vane.  Thus  these  resistance*  are  initially  very  large, 
decrease  rapidly  to  plateau  value*  aa  the  motor  reaches  full 
thru*!,  and  increaae  again  during  burnout.  This  behavior 
has  bean  modeled  by  postulating  that  the  film  coefficients 
begin  and  end  at  IX  of  the  full-flow  values  and  follow  the 
ramp-up,  plateau,  raap-dcvn  profile  of  the  thruat  schedule. 

For  comparison  with  the  simulation  reaulta,  the  vene 
temperature  data  were  likewise  converted  to  a  continuous 
record  by  means  of  the  transfer  function  furnished  by  NWC 
(131.  Thia  transfer  function,  which  was  deduced  from  the 
teat  data,  closely  approaches  the  actual  measurements  and  is 
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given  as 


THRUST  <1.  233s*l> <50. 76s  +  l> 


<  In  this  expression,  thrust  in  Newtons  gives  a  temperature 
response  in  degrees  Fahrenheit  above  ambient. ) 

Vane  Thermal  Response 

Figure  4  shows  the  computed  vane  temperature  history 
using  "current  best  estimates '  o£  the  problem  parameters 
described  previously.  This  fi-:  re  shows  that  the  predicted 
temperature- time  history  r  the  reference  shaft  node  is  well 
in  excess  of  that  obtained  in  the  tests. 

In  the  model  a-  it  is  presently  configured,  the  quanti¬ 
ties  encumbered  by  the  greatest  amount  of  uncertainty  are 
those  associated  with  the  boundaries  of  the  thermal  system 
--  the  convective  inputs  and  the  conductive  cooling  of  the 
van©  by  the  thermal  sink  effect  of  the  mount.  Some  prelimi¬ 
nary  experimentation  with  the  simulation  has  led  to  the 
insight  that  th;?  chief  control  over  the  maximum  shaft 
temperature  is  the  mount  thermal  resistance.  On  the  other 
hand,  the  dynamic  temperature  rise  is  mainly  affeated  by  the 
aonveative  thermal  resistances  acting  at  the  vsne  surface. 
These  observations  lead  to  the  premise  that  the  effective 
time  constants  far  these  two  processes  are  widely  separated 
--  a  clue  to  which  is  given  in  the  previous  transfer-func¬ 
tion  model. 

Adjustment  of  these  factor.’,  leads  results  of  the  sort 
shown  in  Fig.  5.  which  illustrates  the  sensitivity  of  the 


shaft  temperature  response  to  the  convective  film  coeffi¬ 
cients.  It  will  be  seen  from  this  figure  that  in  order  to 
obtain  reasonable  agreement  between  theory  and  experiment  it 
is  necessary  to  reduce  these  coefficients  (with  correspond¬ 
ing  increases  in  the  respective  thermal  resistances)  by  a 
factor  on  the  order  of  6054  to  80)4  --  the  ratio  FR  in  Fig.  5. 
The  need  for  this  offset  is  largely  attributable  to  omission 
of  the  cooling  effects  of  radiation  and  ablation.  It  is 
apparent,  however,  that  the  model  is  capable  of  reproducing 
the  main  transient  features  (time  constants)  of  the  vane 
thermal  response. 

Figure  6  illustrates  the  sort  of  agreement  with  test 
dat*  that  is  achievable  using  the  present  model.  To  achieve 
the  results  shown  in  Figs.  5  and  6, result,  a  mount  thermal 
resistance  of  R3Q  »  10.3  K/W  wos  uaed  --  a  value  not  incon¬ 
sistent  with  the  mass  and  configuration  of  the  mount.  Whe¬ 
ther  or  not  the  reduction  factor  of  72%  is  consistent  with 
the  effect  of  radiative  and  ablative  oooling  is  yet  to  be 
determined.  In  any  case,  it  ie  noteworthy  that  the  agree¬ 
ment  illustrated  in  Fig.  6  ia  attainable  by  means  of  a  sin¬ 
gle  constant  factor. 

In  addition  to  providing  the  framework  of  a  working 
aosputational  code,  the  results  of  the  simulation  indicate 
that  it  may  be  feasible  to  predict  the  bulk  thermal  behavior 
of  selected  aritical  elements  of  the  Jet  vane.  In  order  to 
obtain  these  preliminary  results,  no  particular  attention 
has  been  paid  to  the  precision  with  which  the  boundary  ther- 


Figure  6.  Temperature  vs.  time 


mal  resistances  are  calculated.  For  the  purposes  of  this 
work,  the.  aitnplicity  of  the  model  gives  some  encouragement 
that  it  may  be  adaptable  to  the  deductive  or  system-identi¬ 
fication  mode  of  analysis.  Preliminary  results,  reported  in 
subsequent  sections  of  this  report,  indicate  that  this  is  in 
deed  the  case. 

VI.  STATE  SPACE  MODEL  FORMULATION 

The  results  presented  thusfar  have  been  obtained  from  a 
sequential  integration  of  the  energy  balance  relationships 
at  each  calculation  node,  and  Eq.  <  2 )  illustrates  the  form 
ip  which  these  relationships  have  been  oast  in  the  DSL  pro¬ 
gram  for  purposes  of  simulation.  If  the  physical  quantities 
< resistances  and  capacitances)  appearing  in  the  governing 
system  of  equations  are  truly  constant,  then  the  system  is 
linear  and  tha  »f f ioiency  of  a  state  space  formulation  may 
be  brought  to  bear.  The  matrix  formats  that  characterize 
the  state  spaoto  formulation  are  particularly  useful  in  pro¬ 
viding  a  systematic  presentation  of  the  mathematical  model. 

In  addition,  the  following  results  have  been  obtained 
using  the  personal  computer  version  of  the  software  package 
called  •HATHl&x*  (product  of  Integrated  Systems  Incorpor¬ 
ated,  Palo  Alto,  Calif.).  Incorporated  in  MATRIXx  is  a  vast 
variety  of  matrix  manipulation  schemes  that  are  often  found 
to  be  essential  in  the  modeling  and  simulation  of  dynamic 
systsms  (particularly  control  systems).  Included  are  means 
for  conversion  from  continuous  to  discrete  systems  and  for 
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transformation  of  systems  from  state  space  to  transfer  func¬ 
tion  models  and  vice-versa.  The  system  identification  capa¬ 
bilities  of  MATRIXx,  essential  to  the  goals  of  this  study, 

\ : ... 

4  are  discussed  in  a  subsequent  chapter. 

Mathematical  relationships  are  expressed  in  state  space 
form  by  arranging  them  such  that  the  highest  order  deriva¬ 
tives  of  the  dependent  variables  (the  temperatures  at  the 
various  nodes,  in  this  case)  are  given  as  linear  functions 
of  the  lover  order  derivatives.  For  node  2,  for  instance. 

Eg.  (1)  may  be  written 


T2  * 


la  *  _i_  1 

C2  RF2  Ri2  R23 
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Referring  to  Fig.  3,  the  corresponding  relationships  for 
nodes  1  and  3  are  obtained: 
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The  internal  node  designated  (i)  in  Fig.  3,  which  was 
included  to  allow  separate  estimation  of  the  tip  and  vane 
resistances,  is  not  a  storage  node  and,  aa  previously  noted, 
energy  balance  for  this  node  doea  not  Include  a  rate  term: 

Ti  -  Ti  Ti  '  T2 

R a  -  <  6 ) 

Rii  Ri2 

Equation  <6>  aay  be  aolved  for  and  Eqs.  <3)  and  (4)  may 
be  used  to  obtain  an  expression  for  its  rate  of  chu.ige  in 
the  format  of  Eqs.  <3>-( 3).  This  would  lead  to  a  four-state 
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1  formulation  for  the  four  node  temperatures.  For  efficiency 
V  Of  computation,  however,  it  seems  good  practice  to  eliminate 
such  non-storage  nodes  from  consideration  --  temperatures  at 
such  nodes  can  always  be  calculated  from  algebraic  expres- 
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sions,  such  as  Eq.  (6),  involving  the  temperatures  at  adja¬ 
cent  nodes.  Using  Eq.  <6>  to  eliminate  from  Eqs.  (3)  and 
(4)  gives  the  following  3-state  sett 
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Arranged  in  such  a  fashion,  it  is  apparent  that  the 
nodal  expressions  possess  s  certain  amount  of  symmetry. 


Each  of  the  coefficients  of  tesperature  on  the  righthand 
sides  of  Eqs.  <7>-<9>  have  the  dimension  of  inverse  time 
and,  in  fact,  the  RC  products  are  representations  of  the 


tins  constants  describing  the  energy  transport  proci 


occur ing  at  and  around  the  nodes.  To  further  illustrate 


this  formal  srrcngeaent.,  it  is  useful  to  define  the 
following  parameters  (inverse  tise  constants): 
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bu=  i/ClRF1 


b22  =  1/C2RF2  ' 


b33  =  1/C3R3G 


With  thfSS  definitions,  Eq3.  <7>-(9>  may  be  written 

i=  AT  ♦  |U  (11) 

where  the  "state*  consists  of  the  n  a  3  nodal  temperatures: 

I  =  <T1  T2  Ta>'  , 

the  "input*  vector  contains  the  m  *  3  boundary  temperatures: 
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and  tha  nxn  matrix  A  and  the  nxra  matrix  B  are  given  by 


4 


Equation  < 11 )  is  expressed  in  the  state-space  form  familiar 
to  the  analysis  of  controls  systems,  and  if  the  vector  of 
input  temperatures  y,  are  specified  the  time  histories  of  the 
nodal  temperatures  are  easily  obtsined. 

Because  the  "ground*  temperature  is  a  specified 
constant  in  this  nodal  (in  tha  case  at  hand  it  is  290K)  the 
calculation  may  be  simplified  by  referring  all  temperatures 
to  ground.  This  has  the  effect  of  removing  the  ground  tem¬ 
perature  as  an  input  (it  is  zero  relative  to  itself)  and  the 
input  vsctor  becomes  y.  ■  lTri  T^)*,  With  only  the  recov¬ 
ery  temperature*  as  inputs  (m  *  2)  the  laat  column  in  the 
input  coefficient  matrix  §.  ia  redundant  and  must  be  deleted. 

The  three  nodal  temperatures  are  "outputs, ■  in  a  con¬ 
trol  systems  sense,  and  the  usual  output  relationship  may  be 
written  as 


where  C  is  a  3x3  identity  matrix  and  D  is  a  3x2  null  matrix. 
As  a  final  step  prior  to  simulation,  the  four  matrices  may 
be  combined  into  a  single  partitioned  matrix,  S,  containing 
the  four  submatrices  A,  B,  C,  and  D.  This  is  particularly 
useful  when  dealing  with  the  MATRIXx  software.  Thus 
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In  this  formulation  the  thermal  characteristics  of  the 
eyetea  are specified  in  the  4  matrix  and  the  input  resis¬ 
tances  -are  found  in  the  elements  asking  up  the  &  matrix. 

Note  that  the  diagonal  elements  of  the  4  matrix  are  the 
negative  inverses  of  the  characteristic:  time  constants  for 
the  thermal  response  of  eech  of  the  nodee.  Ae  defined  in 
Eq«  <10 > ,  these  diagonal  element a  are  sums  of  the  values  on 
the  corresponding  row,  including  the  input  matrix.  This 
relationship  expresses  the  conservation  of  energy  embodied 
in  the  analysis  and  the  diagonal  elements  cannot  be  indepen¬ 
dently  varied  without  violating  this  principle.  This  is  an 
important  conatraint  when  ayatea  identification  procedures 
are  applied  to  the  system. 


Th©  generality  of  the  formulation  ia  apparent,  and 
there  would  be  no  particular  difficulty  in  adding  storage 
nodes  to  the  model  in  order  to  obtain  a  more-detailed  view 
of  the  temperature  distribution  within  the  vane.  Such 
refinement  at  this  point  seems  unwarranted,  however,  since 
little  information  is  available  to  allow  a  concomittant 
refinement  of  the  distribution  of  the  thermal  inputs.  In 
fact,  the  present  model  would  not  be  drastically  affected  if 
the  stagnation  point  and  boundary  layer  heat  transfer  were 
combined  into  a  single  input,  thus  reducing  the  model  to  a 
tvo-otate  system. 


SIMULATION  IN  STATE  SPACE 

Figure  3  includes  a  table  of  the  thermal  resistances 
and  nodal  capacitances  that  have  been  estimated  for  the 
problem  at  hand.  (In  estimating  the  "effective*  input 
resistance*  given  in  Fig,  3,  the  calculated  heat  transfer 
coefficients?  have  been  reduced  by  the  721  factor  leading  to 
the  results  described  in  the  previous  chspter.  >  When  these 
values  sre  incorporated  into  the  present  snalysis,  th#  A  and 
8,  matrices  take  the  following  values? 
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0.  18  8  - 
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For  the  tip  node  the  input  temperature  ia  the  stagna¬ 
tion  temperature  (2650K)  and  for  th#  vane  the  recovery  tem¬ 
perature  calculated  by  the  methods  of  Ch.  IV  ia  2S50K  (very 
nearly  the  same  aa  the  stagnation  temperature  in  this  high- 


speed  compressible  flow).  When  referred  to  an  ambient  tem¬ 
perature  of  290K,  the  appropriate  input  values  for  tempera¬ 
ture  are 

Tr1  =>  2360K  and  TR2  =  22S0K 
Although  the  thruat  waa  used  aa  the  fundamental  input  for 
the  simulations  previously  described,  these  temperatures  may 
be  assumed  to  follow  an  input  schedule  that  is  mare-or-leaa 
synchronous  with  the  thruat.  Accordingly,  in  the  simula¬ 
tions  that  follow  the  recovery  temperatures  have  been  given 
the  same  ramp-up,  plateau  (at  the  values  given  above),  and 
ramp-down  contour  that  is  described  in  Ch.  V. 

Lineir— T.iB!B-inviur.lmnt_(.L.TI-)-Mo.del 

With  the  model  matrices  values  given  in  Eq.  (13)  and 
the  scheduled  input  temperatures  described  above,  simulation 
of  the  vane  thermal  response  is  accomplished  using  the  HAT- 
KIXx  L5IN  command  C14,  p. 10-63.  In  thia  case,  the  syntax  is 

[T,XJ  -  LSIff  ( g,  NS,  y,  DELTAT > 

where  §.  ie  the  systes  matrix  formed  as  in  Eq.  (12)  and  in¬ 
corporating  the  values  given  in  Eq.  (13);  NS  is  the  number 
of  states  (NS  «  3  --  the  three  nodal  temperatures);  U.  is  the 
input  err ay,  in  thia  case  expressing  the  scheduled  input 
tespereturss )  and  DELTAT  is  the  specified  time  increment 
between  steps  in  the  j£  array.  In  all  simulation  cases  de¬ 
scribed  herein,  a  tiae  increment  of  0. la  haa  been  used  for  a 
duration  of  15e.  Thus  the  time  vector  haa  dimensions  151x1 
and  the  y.  array  is  151x2  for  the  two  inputs. 

Figure  7  showa  the  reaulte  of  the  elaulatlon  assuming 


thermal  response  —  LTI  baseline  simulation 


that  the  elements  of  the  A.  and  B  matrices  are  constant  with 
time.  Aa  expected,  the  results  are  unsatisfactory  because 
the  linear  'time-invariant  (LTI)  assumption  does  not  account 
for  the  variation  of  the  heat  transfer  coefficients  with 
time.  In  other  wards,  the  elements  of  the  B  matrix  and  the 
diagonal  elements  of  the  A  matrix  must  be  calculated  at  each 
time  step  in  the  simulation. 

Linear  Time- Variant  (LTV)  Model 

With  HATRIXx,  the  accomodation  of  time-varying  parame¬ 
ters  may  be  relatively  easily  aahieved  by  using  the  SYSTEM 
BUILD  capability.  SYSTEM  BUILD  provides  an  interactive, 
menu-driven  graphical  environmsnt  lor  building,  modifying, 
and  editing  computer  simulation  models.  Any  combination  of 
linear,  non-linear,  continuous- time,  discrete-time,  or  mul¬ 
ti-rate  models  that  describe  a  system  may  be  constructed 
from  a  library  of  basic  building  blocks.  Systems  may  be 
modeled  by  dividing  the*  into  individual  components  with 
each  component  being  described  by  a  specific  type  of  func¬ 
tional  block.  Super-blocks  can  be  ueed  to  represent  as¬ 
semblies  af  individual  blocka  and  the  helrarchioal  design 
provides  an  organised  representation  of  the  physical  system 
that  is  Modeled, 

In  SYSTEM  BUILD,  any  super -block  or  set  of  nested  su¬ 
per-blocks  can  be  aieuXated,  linearized,  and  analyzed.  Ho 
user  coding  is  necessary  (although  it  is  allowed),  and  SYS¬ 
TEM  BUILD  can  be  driven  by  uaer -defined  coaaand  procedures 
that  autoaatically  build  or  modify  system  models.  All  of 


this  capability  is  available  for  PC  hardware  (which  has  been 
used  thusfar  in  this  study)  and  the  net  effect  of  the  SYSTEM 
BUILD  capability  is  that  the  user  may  concentrate  efforts  on 
design,  analysis,  and  simulation  in  an  efficient  graphical 
environment. 

The  SYSTEM  BUILD  Model.  In  this  study,  the  vane  ther¬ 
mal  behavior  has  been  modeled  by  means  of  seven  super-blocks 
named  as  follows: 

N0B1IN,  NQD2IN,  N0D3IN 
N0DE1,  N0DE2,  NQDE3 
VANES 

The  appendix  provides  illustrations  of  these  super-blocks, 
and  a  brief  discussion  will  be  provided  here. 

The  first  three  blocks  are  compute  the  parameters  ne¬ 
cessary  for  the  state  equation.  The  input  to  both  NODI IN 
and  NQD2IN  is  the  ramp-up,  plateau,  ramp-down  schedule  pre¬ 
viously  described,  but  vith  a  plateau  value  of  unity.  In 
NQD1IN,  for  example,  this  input  is  multiplied  by  TR1  and 
as  gains  to  form  the  time-varying  element  in  the  first  row 
of  the  input  coefficient  array  B.,  Also  in  MODI  IN,  the  quan¬ 
tity  0^2  generated  by  means  of  the  step  function  and 
this,  in  turn,  is  used  to  form  -a.  ■  a,  +  b.  ,  aa  given  in 

Eq.  (10),  The  outputs  of  NCD1IN  are  -a^  ,  bliTRl'  ®nd  ai2' 
N0D2IN  and  N0D3IN  are  similar  in  function,  with  the  excep¬ 
tion  that  N0D3IN  does  not  require  the  time-varying  input. 

The  first  three  super -blocks  listed  above  ere  elements 
of  the  second  three,  respectively.  In  NQDE1,  for  example. 


the  external  input  is  combined  with  the  outputs  of  NODI IN 

to  form  the  quantity  a11T1  +  ai2T2  *  bllTRl  =  ^1‘  An  inte* 
gration  is  also  within  N0DE1  to  obtain  from  f  .  T  is 
the  output  of  the  N0DE1  3uper-block  and,  in  a  similar  way, 
N0DE2  and  N0DE3  provide  outputs  and  T^. 

Super-block  VANE5  •wires"  N0DE1,  N0DE2,  and  N0DE3  to¬ 
gether  in  order  to  provide  simultaneous  solution  for  the 
three  nodal  temperatures.  In  VANES  (see  appendix),  N0DE1 
receives  the  dual  ramp  input  and  while  generating  T, > 
N00E2  receives  the  dual  ramp  input,  T^,  and  while  gener¬ 
ating  T2*  and  N0DE3  receives  and  outputs  T3> 

To  execute  the  SYSTEM  BUILD  model,  the  program  is  first 
checked  by  means  of  the  ANALYZE  command  Cl5],  Following 
this,  one  of  six  integration  methods  may  be  selected  from 
the  menu  provided.  The  results  described  here  have  been 
obtained  with  a  variable-step  Kutta-Merson  scheme  which  is 
an  explicit  fourth -order  method  that  employs  the  largest 
tide  step  possible  while  remaining  within  error  tolerances 
(maximum  etep  size  is  the  time  inarement  specified  for  the 
simulation). 

Simulation  is  achieved  by  executing  the  SYSTEM  BUILD 
commend  SIM,  as  in  *  SIH(T,y.>",  where  the  time  vector  T 
and  the  input  y.  must  be  provided.  In  the  case  at  hand,  T  la 

the  151x1  vector  containing  IS  seconds  of  0. 1-second  steps 
and  U  is  the  151x2  array  exhibiting  the  raop-up,  unity  pla¬ 
teau,  ramp  down  schedule  previously  described. 

Figure  8  shows  the  results  of  simulation  ualng  the 


Figure  8.  Vane  thermal  response  —  LTV  baseline  simulation 


: SYSTEM  BUILD  model  to  provide  linear  time-variant  (LTV) 
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;£w$P'.'V  behavior  of  the  vane  heat  transfer  coefficients.  In  this 

x  v  figure  the  improvements  attained  with  the  SYSTEM  BUILD  LTV 

iXy  model  are  apparent.  With  the  model  now  cast  in  this  format 

it  is  possible  to  take  the  system  identification  approach 
that  has  been  the  ultimate  goal  of  this  study. 


VIII.  PARAMETRIC  SYSTEM  IDEHTIFICATIOH 

The  philosophy  of  parameteric  system  identification 
<PSX>  is  to  use  measured  response  data  to  "massage*  the 
parameters  of  a  system  model  such  that  simulation  results 
match  those  of  experiments.  Some  of  this  hss  already  been 
accomplished  by  the  menual  variation  of  parameters  described 
An  Ch.  V,  but  PSI  will  now  be  used  to  complete  this  process 
in  s  systematic  way. 

In  the  problem  at  hand,  it  ia  reckoned  that  the  most 
uncertain  parameter  values  ere  those  associated  with  the 
heat  transfer  processes  et  the  vane  surface  and  the  thermal 
sink  effects  of  the  vane  mount.  The  vane  temperature  to  be 
matched  is  that  of  the  shaft  <T^>,  and  although  the  stagna¬ 
tion  heat  transfer  process  at  the  vane  tip  has  been  modeled 
this  provides  relatively  little  heat  flux  into  the  vane  and 
is  therefore  of  secondary  importance  as  far  as  the  thermal 
response  at  ths  shaft  is  concsrned. 

Thus  the  results  dsscribed  hsrs  have  bean  obtained  by  a 
two-parameter  variation  involving  b02  and  a^.  It  is  to  be 
recalled  that  the  variation  of  will  also  require  a 


change  in  a^  in  accordance  with  the  energy  balance,  and 
that  this  is  accomodated  in  the  SYSTEM  BUILD  model.  In  ad¬ 
dition,  the  value  of  has  within  it  which  is  not  used 

in  the  simulation  since  the  relative  ground  temperature  is 
zero.  Nevertheless,  a  knowledge  of  a  will  lead  to  the 

w  *-J 

appropriate  value  of  b^  since  b^  =  -a,^  “  a32  a8  given  in 
Eq.  (10). 

The  MATRIXx/PC  family  includes  a  variety  of  capabili¬ 
ties  that  provide  for  data  analysis,  parametric  and  non-par- 
metric  system  identification,  and  digital  filtering.  In 
this  work,  PSI  is  achieved  by  using  a  scheme  known  as  "maxi¬ 
mum  likelihood  estimation"  (MLE)  C 1 6  3 .  As  its  name  implies, 
MLE  maximizes  the  likelihood  of  the  parameter  estimates 
given  the  observations  (in  this  case,  the  measured  transient 
response  of  Tg>.»  MLE  operates  on  an  entire  time  history 
(rather  than  one  sample  at  a  time)  and  is  applicable  to 


multiple'’ input /multiple -output  (MXMQ)  systems  auah  as  that 
uadpr  analysis  here. 


.  .  The  command  that  executes  the  MATRXXx  version  of  the 
MLE  procedure  is  NAXLIKE  C173,  with  the  following  syntax: 

£ YHAT,  PI  »  MAXLIKE  <  |J,  T3,  P0,  *  VBM2  * , NIT  > 

Items  to  be  provided  by  the  user  ere  U,  containing  the  time 
vector  and  array  of  input  coefficients}  T3,  the  vector  of 
observed  values  that  are  the  standards  for  comparison;  PQ, 
the  vector  of  initial  estimates  for  the  parameters  to  be 
varied ; 'VBM2**  a  user-supplied  command  file  that  creates  the 


values  of  T3  based  upon  U  and  the  current  values  of  P;  and 
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NIT,  the  specified  number  of  iterations  to  be  performed. 

The  results  of  MAXLIKE  are  YHAT,  the  system  response 
.Vising  current  values  of  the  parameters;  P,  which  contains 
th®»®  current  values?  and  RSS  which  is  the  sum  of  the 
.squares  of  the  differences  between  the  measured  and  pre- 

dieted  values  at .each  time  step  --  RSS  =  sum(T3  -  YHAT)  . 

£•  ’ 

Aa  has  been  shown,  the  LTV  nature  of  the  vane  heat 
transfer  problem  requires  the  use  of  SYSTEM  BUILD.  Thus,  in 
-the  problem  at  hand,  the  command  file  'VBM2'  reads  as  fol- 

$T.. 

lovai 

BUILD,  EDIT,  N0D2IN,  0,  EXAMINE,  2,  NEW,  P ( 1 ),  DOME,  TOP,  ..  . 

EDIT, N0D3IN, 0, EXAMINE,  3,  NEW,  P<2),  DONE,  TOP,  .  .  . 
ANALYZE,  VANES 
Y3*SIM<U<  j,  1>,U(  i,2:3)  >  j 
Y;HAT»Y3<  t,3)? 

RETURN 
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i,  %  •  Although  this  language  may  not  be  familiar  to  the  reader,  it 
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.?  may  be  poeaible  to  deduce  that  in  the  SYSTEM  BUILD  lexicon 
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these  inatruotiona  sequentially  update  P(l)  (*  b22 )  which  ia 
in  poeition  2  of  super-block  N0D2IN  and  P(2)  < *  a33>  which 
la  in  position  3  of  super-block  NQD3IN.  Y3  ia  the  name  gi¬ 
ven  to  ths  three-temperature  array  obtained  from  simulation 
using  the  current  parameter  values,  P(l)  and  p<2>,  with  YHAT 
(«  T^>  contained  in  third  column.  It  should  be  noted  that 
in  each  iteration  the  MAXLIKE  procedure  executes  ’ VBM2 *  at 
once  for  each  parameter  value  under  consideration. 

Figure  9  is  the  result  of  a  PSX  run  of  three  iterations 


Figure  9.  Vane  thermal  response  —  LTV  simulation  after 
system  identification. 


using  initial  parameter  values  given  previously  <fc>22  =  0.20 
and  a^g  3  ~0. 48).  The  "maximum  likelihood"  values  obtained 
through  the  use  of  PSI,  and  giving  the  results  shown  in  Fig. 
9,  are  2027  and  a33  =  -0. 4565.  The  value  of  the 

residual  associated  with  these  results  is  RSS  =  59.4.  (This 
is  a  very  small  value,  bearing  in  mind  that  it  is  a  sum  of 
151  squared  temperature  differences. ) 

Diacusalon  of  the  PSI  Results 

As  vu  done  in  Ch.  V,  the  experimental  shaft  tempera- 
ture  data  (shown  as  x's  in  Pigs.  7-9)  have  been  obtained 
from  the  transfer  function  provided  by  NWC.  In  this  case, 
the  temperature  above  ambient,  in  Kelvins,  is  given  by 


4.  48 


(14) 


THRUST  <i.233s-l)<50.76a-l) 
where,  as  before,  the  thrust  has  the  dimension  of  Newtons 
(see  Fig.  6).  h  state-space  equivalent  of  this  expression, 
which  is  easily  obtained  by  means  of  the  MATRIX*  function 
SFORN,  is 

k  m  *  B.u 

.  y  *  •  Ru 


where 

A  - 


'  -0.8307  -.0.0160 

'  1  ' 

*  &  a 

.  i  0  , 

.0  , 

g  *  to  0.0716) 


and  £  *  0.  <Xr»  this  representation,  y  ia  the  shaft  tempera¬ 
ture  Tg  end  u  is  the  thrust.)  In  the  combined  matrix  format 
previously  discussed  in  connection  with  Eq.  (12),  these  ma¬ 
trices  yield 


s  * 


-0.6307  -0.0160  1 

1  0  0 
,  0  0. 0716  0  . 


<  15) 


The  agreement  illustrated  in  Fig.  9,  which  is  clearly 
excellent,  can  be  evaluated  in  several  ways.  For  instance, 
the  temperature  response  resulting  from  the  identified  sys¬ 
tem  may  be  cast  in  a  form  that  is  directly  comparable  with 
that  provided  by  NWC.  To  accomplish  this,  th©  system  model 
may  be  converted  to  the  2-state  matrix  form  given  above,  and 
the  three  parameters  of  S  subjected  to  system  identifica¬ 
tion,  This  has  been  done  by  means  of  the  MAXLIKE  procedure 
with  the  following  result: 


0.  6216  -0.0173  1  * 

1  0  0 
0  0.0716  0 


with  a  value  of  RS8  *  63.4.  Comparison  of  this  result  with 
Eg,  <15)  demonstrate*  the  agreement.  Using  the  TFQftM  com¬ 
mand  of  HATRXXx,  the  atate-apac®  representation  may  be  con¬ 
verted  back  to  tranafer-funotion  form,  yielding 


T  * 

*3  4.  13 _ __ 

THRUST  *  < 1. 2S0a*l) ( 46. OQa^l ) 


Comparing  with  Eq.  <X4>  the  agreement  ia,  again,  remarkable. 
The  reader  should  bear  in  mind  that  these  arguments  are  made 
here  only  for  the  purpose  of  demonstrating  the  effectiveness 
Of  the  P€)X  procedure  in  this  case.  The  various  transfer 
functions  and  parameters  in  these  models  have  no  relation  to 
the  physical  situation  --  in  feat,  since  they  represent  LTI 
systems,  they  are  wholly  Inadequate  for  any  purpose  other 
then  expressing  the  experimental  measurements  simulation 


results  In  functional  forms. 


Several  sets  of  initial  values  have  been  used  in  iden¬ 
tifying  the  parameter  values  in  the  physical  model,  and  the 
results  have  been  found  to  be  independent  of  this  factor. 

In  addition,  a  number  of  teats  have  been  performed  to  deter¬ 
mine  the  effect  of  varying  other  parameters  in  the  system 
model  and  no  signif icsmt  improvement  hae  been  obtained  on 
the  results  illustrated  in  Fig.  9.  The  selection  of  the 
appropriate  parameters  to  vary  in  the  PSI  process  requires 
some  appreciation  for  the  phyaical  process  modeled.  Thus  it 
has  been  observed,  for  instance,  that  the  selection  of  too 
many  parameters  sill  lead  to  unrealistic  values,  such  as 
negative  thermal  conductances,  unless  constraints  are  added 
to  the  HA XL IKE  procedure.  These  and  other  aspects  of  the 
method  are  under  continued  study. 

XX,  CONCLUSIONS  AND  RECOMMENDATIONS 
The  eaphaais  in  the  work  reported  here  has  been  to  in¬ 
vestigate  the  feasibility  of  using  system  identification 
methods  to  develop  tools  for  use  in  the  design  of  jet  vane 
TVC  devices.  Before  proceeding  to  conclusions  in  this 
regard*  ho* ever,  it  may  be  useful  to  comment  on  the  results 
obtained  in  the  case  that  has  been  analyzed. 

The  excellent  agreement  between  measured  and  predicted 
Shaft  temperature  reaponaea  shown  in  Fig.  9  gives  some  cre¬ 
dibility  to  the  model  parameters  used  in  the  simulation. 


tu 


Thessf  may  In  turn  be  related  to  vane  physical  properties  to 
lend  some  insights  regarding  those  quantities  that  were 
initially  subject  to  relatively  large  uncertainties.  For 
the  heat  sink  effect  of  the  shaft  mount,  for  instance,  the 
parameter  b^  =  ”  “  a32  ma^  °tt.'Vined  from  the  identi¬ 

fied  value  of  a^g  -  -0.46  and  the  assumed  (and  relatively 
certain)  value  of  a^  =  0.38.  The  result  is  b33  =  O.Ofi  s  1 
and  this  may  in  turn  be  related  to  the  van©  physical  proper¬ 
ties  since  b__  =  1/C,.,R__.  This  new  value  is  down  from  the 
33  3  3G 

original  value  of  0.  10  s  *  (Fig.  3),  and,  since  RC  -  mcL/kA, 
the  result  may  be  used  to  adjust  the  rationale  used  in  esti¬ 
mating  the  effective  mass  and  geometry  of  the  thermal  sink 
representing  the  mount  (assuming  that  the  heat  capacity  and 
thermal  conductivity  are  firm  values). 

The  other  parameter  identified  in  the  FSI  analysis  is 


the  vane  thermal  input  coefficient,  »  0,2027  s 


-l 


If 


this  value -ie  used  to  recalculate  the  convective  heat  trans- 

2 

fer  coefficient*  a  value  of  h  «  638  W/m  K  is  obtained  which 
is  about  72*  lees  than  that  eetimated  in  Ch.  IV  for  turbu¬ 
lent  boundary  layer  Convection.  The  conclusion,  ae  has  been 
mentioned#  is  that  the  complex  processes  of  ablation  and 
rsdiatlda  account  for  a  considerable  amount  of  cooling  of 
the  vane.  It  is  Important  to  keep  in  mind  that  the  single 
values  quoted  here  are  useful  only  inasmuch  as  they  give 
good  results  for  predicting  shaft  temperature  response  using 
the  model  that  has  been  constructed.  In  this  case,  the  mo¬ 
del  includes  a  ramp-up,  plateau*  ramp-dovn  behavior  for  b^ 


so  that  the  above  result,  which  applies  to  the  plateau  value 
of  the  "effective"  heat  transfer  coefficient,  accounts  in  a 
crude  way  for  the  fact  that  vane  heat  transfer  effects  vary 
widely  throughout  the  firing.  The  assumption  of  a  constant 
plateau  value  for  h  is,  of  course,  subject  to  question  and 
refinement  if  sufficient  information  is  available. 

Referring  to  the  vane  surface  temperature  responses 
shown  in  Fig.  9,  it  is  seen  that  a  maximum  of  about  1240K 
OSSGK  above  ambient)  is  reached  at  about  3.5  seconds  (3  se¬ 
conds  after  motor  ignition).  This  value  is  well  below  the 
melting  ooint  of  the  vane  material,  and  the  tensile  modulus 
is  still  about  22,000  kai  at  this  temperature  (though  it  is 
w« 11  below  the  ambient  value  of  about  45,000  ksi).  It  might 
bs  predicted  from  these  results  that  the  1/ 4-scale  vane  that 
has  been  Modeled  would  remain  more-or-less  intact  during  a 
firing  of  the  type  that  has  been  simulated.  A  noteworthy 
feature  of  these  developments  is  that  it  would  be  a  rela¬ 
tively  simple  matter  to  simulate  the  thermal  response  of  a 
full-scale  vane.  Each  of  the  paraaetere  in  the  present  mo¬ 
del  is  amenable  to  direct  scaling  procedures  such  as  those 
that  have  been  given  in  Ch.  IV. 


The  results  thuafar  obtained  in  thia  study  justify  some 
confidence  that  the  thermal  behavior  of  the  jet  vane  can  be 
modeled  using  a  ralatlvely  simple  and  etraight-forward  model 
Structure.  But  although  the  structure  of  the  model  is  well 
defined,  some  of  the  values  of  ths  system  parameters  necea- 


sary  to  make  the  model  "work*  are  subject  to  considerable 
uncertainty.  In  this  study  the  power  of  the  system  identi¬ 
fication  method  has  been  brought  to  bear  in  removing  some  of 
this  uncertainty.  In  addition,  it  is  worth  mentioning  that 
the  relatively  new  capabilities  now  available  with  software 
products  like  KATRIXx  have  greatly  facilitated  these  advan¬ 
ces. 

Future  research  in  these  areas  is  warranted  from  both 
the  modeling  and  system  identification  points  of  view.  Us¬ 
ing  the  model  that  has  been  developed  here,  or  perhaps  an 
even  simpler  model  (without  the  tip  node),  it  is  important 
to  determine  if  the  good  results  obtained  in  the  present 
instance  can  be  expected  in  other  cases.  With  this  in  mind, 
the  results  of  other  firing  tests  should  be  predicted  with 
the  model  and  further  system  identifiastion  conducted  if  it 
proves  to  b©  necessary.  The  basic  issue  would  be  to  deter¬ 
mine  the  robustness  (or  lack  thereof)  of  the  present  model. 

Pro®  the  modeling  point  of  view,  a  number  of  inter¬ 
esting  questions  remain  of  a  more  theoretical  nature.  The 
relationship  between  the  complexity  of  the  model  and  the 
adequacy  of  experimental  comparison  data  ±a  of  some  impor¬ 
tance.  In  addition^  more  realism  could  be  introduced  into 
the  present  aodel  in  an  effort  to  determine  the  sensitivity 
of  the  results  to  such  matters.  For  instance,  it  might  \>e 
possible  to  introduce  a  multiple  thermal  input  model  that 
accounts  for  radiation  and  ablation Haa  well  aa  thermal  con¬ 
vection.  Other  elements  Qf  realism  could  be  introduced  by 


allowing  for  the  temperature  dependency  of  several  of  the 
parameters  now  taken  a a  constant.  These  include  the  thermal 
conductivity  of  the  vane  material  and  the  reference  tempera¬ 
ture  used  in  calculating  the  thermal  properties  of  the  motor 
exhaust  gases. 

Another  aspect  of  the  method  that  may  lead  to  important 
insights  is  that  the  initial  and  final  thermal  responses 
seem  to  be  somewhat  uncoupled  under  the  present  circumstan¬ 
ces*  For  .example*  the  L.TI  response  illustrated  in  Fig.  7  is 
quit®  adequate  during  the  initial  phase.  The  implication  is 
that  the  initial  and  final  transients  might  be  used  to  iden¬ 
tify  particular  parameter  values  under  the  assumption  that 
they  are  time-invariant.  If  this  were  true,  the  handling  of 
nonlinear  behaviors,  should  this  becooue  necessary,  would  be 
greatly  facilitated. 
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Diagrams  of  System  Build  Super-Blocks 
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